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A self consistent field theory that describes a part of a contour loop of a random 
Gaussian surface as a trajectory interacting with itself is constructed. The exponent v 
characterizing the end to end distance is obtained by a Flory argument. The result is 
compared with different previous derivations and is found to agree with that of 
Kondev and Henley over most of the range of the roughening exponent of the random 
surface. 



The study of surfaces, their static statistical properties, as well as growth and 
evolution dynamics has been attracting an ever increasing amount of interest over the 
last two decades. The main reason is that in addition to the many important problems 
of real surfaces arising in condensed matter physics, 1 " 6 this kind of problems is 
intimately related to other problems in physics like directed polymers, 2 ' 7 string 
theory, 8 phase transitions in two dimensions 9 etc. In many of these studies the main 
interest lies in obtaining the distribution functional and the roughness exponent £ 
characterizing these self affine surfaces. The purpose of the present letter is to study 
the statistical properties of contour lines of equal height. In particular I will be 
interested in the exponent v characterizing the end to end distance AT? of a part of 
such a contour of length S, AR S v . The exponent v will be discussed in the case 
of random surfaces that are governed by a Gaussian distribution. This is the simplest 
distribution one can think of, yet it may serve as a reasonable approximation of a real 
surface distribution, in particular, when one is interested mainly in large scale 
properties. 3 

The specific problem I consider in the following was discussed in the past using 
two different methods. The approximate multiscale analysis of Isichenko 10 yields 
V = 7/(10- 3£) , while very powerful arguments by Kondev and Henley 11 yield a 

different result, v = 2/(3 - £). 

The purpose of the present article is to construct a one dimensional field theory, 
that ascribes an energy to any given contour much in the way that it is done for 

1 9 

polymers and self repelling polymers. The main difference is that the "Hamiltonian" 
that ascribes an "energy" to any given contour f(t) is self consistent. Namely, its 

structure involves averages that have to be calculated with respect to the same 
Hamiltonian. The exponent v obtained by this method agrees with that of Kondev 

and Henley 11 but only above £ = j/^ . Below that value the exponent I find from the 

self consistent Hamiltonian is just the end to end exponent of self avoiding walks but 
the use of the self consistent Hamiltonian in that regime is doubtful. 

Consider a self affine surface h = h(r) , where h is the height at the point 
r = (x, y) in the plane. The surface is assumed to be random and governed by a 
Gaussian distribution 

P{h}o, exp [-^fd 2 q q 2(1 ^h q h_ q ], (1) 

where y is the high momentum cut-off and h q is the Fourier transform of h(r) . All 

the Fourier component of h(r) with 1>y a WQ taken to be identically zero. This 

implies, with probability 1, that the surface is smooth on scales less or equal to a . In 
fact, if we extend x and y to the complex variables z x and z 2 , respectively, 
h(z 1 ,z 2 ) is an analytic function in both variables. Note also that since h can be any 
physical scalar field, not necessarily a height, a is the only length in the problem . 
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As a result of the smoothness of the surface, the contour lines are smooth on the same 
scale. Consequently, an arc length can be defined and as any other length it will be 
naturally measured in unite of a . My purpose is to study the statistical properties of 
contours of equal height, f = r(s) , in arc length intervals o < s < S , such that for 

s r ,s 2 in that interval r(s r ) ^ r(s 2 ) . In particular, I will be interested in the exponent 
V , that characterizes the end to end distance 

<[Ar(S)] 2 >~S 2 \ (2) 

where the average < > is taken with respect to the distribution (1). The first goal is 

to obtain a statistical weight for a line starting at an arbitrary point and having the 
parametric representation Ar = Af(s) . Note that the exponent v is not expected to 

depend on a (universality) but non-universal quantities like prefactors will depend on 
a . For example, in eq. (2) above the constant of proportionality on the right hand 
side of the equation is a pure number times a 2(l ~ v) . Setting a equal to zero renders 
S 2 and < Ar 2 > infinite with a proportionality constant that is zero. Therefore, one 
can take the limit a — > only for v but not for S or < Ar (S) 2 > . 

An equal height contour can be obtained by following a "particle" that obeys the 
dynamical rule 

df Vh 

— = zx— = zxu , (3) 
ds IV/il 

where z is a unit vector in the z direction (perpendicular to the xy plane) and u is 
a unit vector in the plane perpendicular to the line. The unit vector z x u is tangent to 
the equal height contour passing through r . The fact that u is a unit vector implies 
that s is indeed an arc length as required. The above "dynamical rule" takes the 
"particle" from any given starting point, along the equal height contour passing 
through that point with "speed" V = 1 . The formulation of the problem becomes more 
tractable if we replace eq.(3) by 

r=- = czxVh, (4) 
dt 

where c = [<l Vh l>] -1 . Now, <l r l> = 1 but we do not expect \r 1= 1 to hold for 
all t . Equation (4) above describes a particle moving with speed 
V(r) =1 Whir) I / <l Vh l> . I will be interested in the following in the exponent 

characterizing the dependence of < Ar 2 > on the "elapsed time" T. How is this 
related to the dependence on S described earlier? It is straightforward to show, 

using the relation S = I V[r(t)]dt , that the average S for fixed T is T. For large T the 

fluctuation in S is expected to become vanishingly small compared to its average. 
This follows from the very plausible assumption, that the correlation 
<(V[r(/ 1 )]-l)(V[r(/ 2 )-l) > vanishes for large \t 2 -t l \, because a large "time" 
difference implies a large distance \r(t 2 )-r(t l )\. If this is the case, a choice of T 
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defines actually a sharp S and the exponent characterizing the dependence of < Ar 2 > 
on T is the same characterizing its dependence on S. A proof that the above 
description is valid in our particular case is still lacking, however, and arguments 
suggesting a broad distribution of S' s for given T can be presented. As long as a 
proof of either of the two possibilities does not exist it is important to note that for our 
purposes a much weaker result is needed. Indeed, even if the speed correlation tends 
to a constant, so that the fluctuation in S is of order T,S scales as T and the 
characterizing exponent is the same. (It is easy but somewhat lengthy to show that the 

2 V 

fluctuation, < AS (T) > n is bound from above by cT T, where c is a constant of 
order 1). 

The probability of a given trajectory (I define one end point to be the origin) 
r(T)= f T ?(f) , (5) 

J o 

is given by 

P{r(t)}=<n6(r(t')-czxVh(f(t'))% T {f(t)} >, (6) 

t'=0 

where % T {r(t)} is a characteristic function that is zero if r(t l ) = r(t 2 ) for some 
t l ,t 2 <T and one otherwise. Since % T does not depend on the specific height 
configuration it can be brought out of the average and then can be described in terms 
of an Edwards interaction 12 so that 

P{? (t) } oc< n 8(7 (t ') - cz x Vh(f (t ')) > x exp{- - C C dt^t^p (t, ) - r(t 2 )] } (7) 

t'=0 ^ Jo Jo 

Introducing a two dimensional vector field, k(t)dn& using the representation 
8{x) = — f e lkx dx , the probability P{r(t)} is found to be 

P{f(t)}oc fDk(t) exp[if T dt k(t) • ?(t)] < exp[-ic f T dtk(t) • (z x Vh)] > x 

J Jo Jo 

xexp{--J o |dt 1 dt 2 5[r(t 1 )-r(t 2 )]} 

where Dk(t) denotes functional integration over the vectorial function k(t) defined 
on the interval < t < T . Since h is a Gaussian random field 

< exp[-z'c fdt k(t) ■ (z x Vh)] >= 

Jo 

(9) 

= exp[- - c 2 f \ T dtdt 2 < k(L ) • £ x Vh[r(L )]k(t 2 ) • £ x Vh[r(t 2 )] >] 

2 Jo Jo 

The needed average can be easily performed 
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< k(t, ) • z x Vh[r(t x )]k(t 2 ) • z X V/*[Vr (f 2 )] >= £ <X ) ■ (I x (? 2 ) ■ (z x V 2 ) 
<MP(0]-M?(g]]^=4(0-(fxV)^ 2 )-(zxV)w[r(r l )-r(g] , (10) 



where Vjacts on functions of r(t x ) and V 2 on functions of r(t 2 ) . 

Equations (8)-(10) define a one dimensional field theory in terms of the fields k(t) 
and r(t) governed by the Hamiltonian 



where A = V 2 w . 

In obtaining eq, (11) I have used zxk — > and rotational invariance. 

The above Hamiltonian presents an exact field theoretic description of the 
problem, that should enable a systematic treatment. In the following I use, however, 

a simplification that will enable the k integration. The approximation I use is to 
replace A[r(t 2 ) - r(t r )] by its average over trajectories rit) . For large values of time 
separations, the distance, Ar =1 r(t 2 )- f{t x ) I, scales as the average over trajectories, so 
that inspite of a fluctuating prefactor of order 1, the replacement of Ar by its average 
still yields a A that has a correct scaling and this is what counts in the simple scaling 
(Flory) argument given in the following. (The fact that A has the correct scaling at 
large time separations is a result of the fact that for large arguments A behaves as a 
power law). 

The remaining question whether indeed the long "time" dependence of Ar is 
dominated only by long "time" separations in the Hamiltonian will be discussed later. 
At present I continue replacing A[r(t 2 ) - f(t r )] by its average for all time differences 



The average [A[[r(^) - r(t 2 )]] is a function of \t x -t 2 \ and going over to Fourier 

transform in t facilitates the k integration. The result is a self consistent 
Hamiltonian for r(t),H{r(t)} , with parameters depending on averages, that have to 
be calculated with respect to the same Hamiltonian 



where r((o) is the Fourier transform of r(t) and (j)((0) is the Fourier transform of 
[A[r(^) - r(/ 2 )]] . For small ft), 0(G)) ftT M with 




(ID 



H{7(t)} = -^\d(0 



2c 




(12) 
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j u = [l-(2-2C)v] 



(13) 



Note that the above Hamiltonian is a generalization of the Hamiltonian of the self 
avoiding polymer that is characterized by ji = . 

I use now a Flory type argument to obtain a relation between Ar and T. This is done 
by scaling: co'=a)T r'=r/Ar and consequently r'(co') = f(o))/(TAr) . As a result 

the first integral on the right hand side of eq. (12) scales as (Ar 2 )/r 1+M while the 

second integral scales as T 2 /(Ar) 2 . In order that the two terms be of the same order 

of magnitude we must have (Ar) 4 oc r 3+M or 



v = (3 + J u)/4. (14) 

It is well known that although the use of a Flory type argument is not an exact 
treatment, it yields in many cases exponents that are extremely close to the exact ones 
and in some cases even the exact exponents. 

Indeed, the exponents obtained from the Hamiltonian (12) for the case ^ = (or 
Li < as will be seen later) is known to be exact. 1314 Using eq. (13) one obtains 

v = 2/(3-0. (15) 

This is the result obtained by Kondev and Henley but only naively so. The reason 
is that, the small (0 (/)((») , is proportional to (0~^ with ji given by eq. (13), as long 
as that ji is positive. If there is no special reason to believe that indeed (f)(0)) can 
vanish for small 0), then when eq. (13) yields a negative ji this means that (j)((0) is 
actually a non zero constant at small 0) (or that ji is actually zero). The end to end 

distance exponent of the self avoiding walk is thus recovered, v = % . 

The result of the above calculation is that the relation v = 2/(3-£) holds only as 

long as £ > . For roughness exponents < £ < v = - 

The origin of the above small discrepancy is unclear yet . Its study as well as a 
more accurate treatment of the problem is postponed to future publications. At this 
stage I am in a position where I can only suggest the direction where the solution must 
be found. In the derivation of the self consistent Hamiltonian (eq. (12)), I have 
replaced A[r(t 2 ) - r(t x )] by its average. Clearly, this is valid only for largely —t 2 \. 
The question is, how is the final result affected by the mistake done for small 
\t y —t 2 I? Since the Hamiltonian H(r,k) (eq. 11) is quadratic in the k's, Gaussian 
integration yielding an exact Hamiltonian that depends only on r(t) is possible. The 

explicit form of the Hamiltonian can be written in terms of A -1 , the inverse of the 
kernel A(t r ,t 2 ) = A[r(t 2 ) - r(^)]. It proves useful to split the kernel into two parts: 
Aj = A[t 1 ,t 2 ]x-6[T -\t 2 -t x I] - and A 2 = A[t 2 ,t l ]x-6[\t 2 -t x \-TJ where 7;isascale 
chosen so as to ensure that for \t 2 -f, l> T a , A[r(t 2 ) - r(^)] may be replaced by its 
average. 
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The Fourier transform of the kernel 




(16) 



/•V /». /». 

is thus the sum of two contributions, A x ((O x ,(0 2 ) , and A 2 (fi) 1 ,G) 2 ) , where A,(ft),,ft) 2 ) 
is bounded and A 2 {(o 2 ,(O x ) is diagonal and has the form A 2 (co 1 ,© 2 ) co~^ 8 , for 

small (0 1 . Therefore, for ji > where A 2 has eigenvalues that diverge as ft), tends to 
zero, A 2 provides an adequate description of the long time behaviour (as given by the 
self consistent Hamiltonian). Note, however, that A 2 dominates over A, , provided 



(l > . The use of the self consistent Hamiltonian is thus restricted to the 
regime g > j4 • The final conclusion is therefore, that the Kondev- Henley relation is 



obtained here only for g > % . It is unclear yet, what the present approach yields of 



The fact that the Kondev- Henley result is supposed to be exact at g = 1115 and the 
naive derivation of their result for all 0<g<l above (eq.(15)) pose intriguing 
problems that deserve further exploration. 
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